SESSION III.- FORUM FOR DISCUSSION AND DEBATE 


The second day was opened by A. Barrie Pittock who made a plea 
for including physical understanding of sources of variances and 
the physical processes in the atmosphere. Statistical models 
looking at the data set without knowing what is going on are likely 
to be misleading. He cited the classic example of water levels in 
Lake Victoria which showed two nice 11-year cycles prior to the 
early 1920s that correlated with sunspots, but then showed much 
shorter, small amplitude cycles until the early 1960s. A massive 
rise of more than a meter then took place and levels have dropped 
only slowly since then. 

To illustrate his point that physical insights can make sense 
of climatic series and provide evidence of causal relationships, 
Pittock showed a time series of precipitation in Seattle with an 
apparent anomalous increase in rainfall in the Puget Sound area 
since 1940. He then showed how this apparent anomaly can be 
accounted for meteorologically with the location of high pressure. 
"Thus," Pittock concluded, "we can use one physical time series 
to account for another." 

Pittock pointed out the high variability of atmospheric ozone 
content, the variance of which changes markedly with altitude 
(Fig. 42), and showed the results of a recent analysis which broke 
the total variance in ozone over Aspendale (38S) down into components 
having different time scales and possible causes (Fig. 43). 

VARIABILITY <nbor) 




Figure 42. Variability of Atmospheric Figure 43. Components of Variance 
Ozone Content in Ozone over Aspendale 
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Noting that spatial patterns of variation give clues to the 
physics behind them, Pittock stressed the need to identify regions/ 
stations which should be monitored to understand apparent trends. 

In specific reference to Hill's methodology, Pittock said, "It 
is not just a matter of selecting equal area boxes but worrying 
about where the boxes are." 

Spatial patterns of mean distributions or patterns of change, 
or eigenvector characteristic patterns, or patterns of correlations 
between stations or with circulation parameters can be used. 

In eigenvector analysis, usually 80% to 90% of the total 
variance can be accounted for by the first eight or so patterns. 

So, Pittock suggested identifying patterns which account for the 
variances, then looking for what might cause them. 

Pittock continued, "A few such patterns usually account for 
most of the variance, leading to physical hypotheses concerning 
causal relationships which can be tested." The dominant patterns 
in many climatic variations are standing waves, due to orographic 
effects and land-sea distribution, and patterns related to the 
strength of the Hadley circulation. These mechanisms, which operate 
on ozone, largely account for correlations between stations (Fig. 44) 
and suggest where monitoring stations should be located. Ozone in 
the southern hemisphere is highly correlated with the latitude of 
the high pressure belt (Fig. 45). 


Correlations R, between Spring (ASO) mean total ozone amounts 
between various pairs of Southern Hemisphere stations. N is the 
number of data pairs, is the percentage of the variance accounted 
for by the correlation, and P is the percentage probability that 
the correlation has occurred by chance. 


STATIONS 

N 

R 

R^ 

P 

Hobart & Wellington v Aspendale 

9 

0.85 

72 

< 1 

Macquarie Isle v Aspendale 

10 

0.81 

66 

< 1 

Amundsen-Scott & Byrd 
V Aspendale 

11 

0.77 

59 

< 1 

Brisbane v Aspendale 

14 

0.73 

54 

< 1 

Darwin v Aspendale 

5 

0.20 

4 

1 arge 

Argentine Island v Aspendale 

8 

-0.33 

11 

large 

Darwin v Brisbane 

• 

5 

0.44 

20 

large 


Figure 44. Correlations Between Stations 
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Time series of amplitudes of 
characteristic patterns should be 
monitored and compared with 
correlated circulation indices. 

A breakdown in well-established 
correlations between ozone varia- 
tions and variations in other 
atmospheric parameters or indices 
would suggest the need to inves- 
tigate anthropogenic causes. 

Hypotheses as to anthropo- 
genic causes should be tested by 
correlating indices of hypothesized 
causes, intermediate effects and 
corollaries, as well as effects on 
ozone. Pittock concluded, "If there 
is a change occurring and there is 
not a change occurring in the 
general circulation, then we'd 
get very suspicious." 
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Figure 45. Total Ozone Versus 
Latitude 


London presented information that the largest variance occurs 
with the largest ozone buildup (in winter), not at the largest total 
amount of ozone, to which Tukey added that "in a system with feed- 
back, arguing with lags is hanging over an abyss because to say it 
occurs is also to say the reverse is true." 

Lovill next discussed a paper. Temporal Variability of Total 
Ozone During 1957-75 , written with his Lawrence Livermore 
Laboratory colleagues Thomas J. Sullivan and John A. Korver. 

The paper, as submitted to the proceedings, follows . 


There are 152 stations that have taken total ozone observations. 
The length of record varies from 6,618 days (July 1957-December 
1975) with observations taken at Aspendale, Australia to as few as 
six days at Woomera, Australia. This paper will use only the data 
from 15 of these 152 stations. The stations were selected on the 
basis of longevity of record and their individual standard 
deviation (a). Each of these 15 stations has a minimum record 
data length of 18 years. The standard deviation of ozone values at 
a station is primarily a function of the instrument calibration and 
daily meteorological variability. 
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We have calculated the standard deviation of total ozone 
variations at a subset of 99 stations. These stations are located 
as shown in Figure 46 (a,b) and their a's are indicated in Figure 47 
(a,b) and Table 10. It is readily obvious that the a's increase in 
value from lower to higher latitudes. The standard deviations range 
from as low as 9 m atm-cm at Huancayo, Peru to as high as 
108 m atm-cm at Yakutsk, U.S.S.R. 

It is worthwhile to compare the a's at stations in similar 
latitude bands in order to obtain an estimation of individual station 
meteorological variability and instrumental accuracy. The a's at 
stations in North America compare well with those in Western Europe 
at selected latitude bands. A comparison of the Western European 
and North American data with those of the Japanese stations also 
indicates similar values as a function of latitude. However a's 
at many stations in the Soviet Union do not compare well with the 
data from North America, Japan, and Western Europe. In the southern 
hemisphere there are considerably fewer stations and the a variability 
is large. Two stations do appear to deviate significantly from the 
average for their latitude band: these are Port aux Francais 

(a = 83 m atm-cm) and Dumont d'Urville (a = 85 m atm-cm). 

Next we looked at regional total ozone variations during the 
18-year period by combining the individual station records for 
selected regions (Figs. 48-51). 

When this is done for the two Canadian stations ozone is observed 
to increase irregularly until 1966; thereafter it irregularly de- 
creases. The combined record of the three Japanese stations indi- 
cates an irregular, slow increase of ozone that is continuing until 
the present. The two Australian stations indicate an irregular 
decrease of ozone continuing until the present. The Indian stations 
show a strong increase of ozone until 1964 and thereafter a slower 
increase and since 1970 a steady amount. 

Next we have expanded our coverage using these 15 stations 
until it is global in extent. We will look at two different 
techniques for analyzing these 15 stations, which we think represent 
the best long-term data record available. In Figure 52 we have 
plotted the 18 years of data from the 15 stations such that each 
station contributes equally. These data in Figure 52, which are 
strongly biased toward the Northern Hemisphere (especially Europe) 
indicate an increase of total ozone until 1970 and thereafter 
a decrease. Figure 53 weights the station data in the Northern 
Hemisphere equally with those from the Southern Hemisphere. In this 
figure it is very difficult to determine a trend of any significance. 
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Figure 46a. Northern Hemisphere Station Locations 



Figure 46b. Southern Hemisphere Station Locations 
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Figure 47a. Northern Hemisphere Total Ozone Observatory Stations 

Used to Calculate Standard Deviations 
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Figure 47b. 


Southern Hemisphere Total Ozone Observatory Stations 
Used to Calculate Standard Deviations 
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Table 10. Ninety-nine stations and their standard deviations 


Station 

Number 

Station Name 

a 

O3 

No. of 
Obser- 
vations 

Lat i tude 

Long! tude 

Station 

Number 

Station Name 

a 

O3 

No. of 
Obser- 
vations 

Latitude 

Long i tude 

3 

Alma Alta 

67 

321 

4521 

43. IN 

76. 5E 

42 

Leningrad 

65 

350 

3997 

59. 5N 

30. 2E 

5 

Dikson Island 

103 

364 

1502 

73. 3N 

80. IE 

43 

Lerwi ck 

60 

35*t 

1(513 

60. IN 

I.IE 

7* 

Kagoshima 

30 

289 

5804 

31. 4N 

130. 3E 

44 

Spi tzbergen 

67 

353 

987 

78. IN 

IS.liE 

8* 

Kodaikanal 

18 

257 


10. IN 

77. 3E 

45* 

Messina 

40 

343 

6342 

38. IN 

15. 3E 

9 

Mount Abu 

17 

255 

2678 

24. 4n 

72. 4e 

47 

Naples 

40 

299 

3139 

40.5N 

|1(.2E 

10* 

New Delhi 

21 

272 

5974 

28. 4n 

77. IE 

48* 

Oxford 

53 

356 

5328 

51. 5N 

l.iw 

n 

Quetta 

30 

277 

3325 

30. IN 

66. 5E 

50 

Potsdam 

50 

347 

3024 

52. 2N 

13. OE 

12* 

Sapporo 

61 

368 

5918 

43. ON 

141. 2E 

5) 

Reykj av i k 

61 

339 

3179 

64. IN 

21. 5W 

13 

Spinagar 

26 

292 

4408 

34. ON 

74.se 

52 

T romso 

75 

337 

2474 

69 . 4n 

18. 5E 

«<.* 

Tateno 

42 

323 

6515 

36. ON 

140. OE 

53 

Uccle 

51 

351 

1411 

50. 5N 

4.2E 

15 

Torishima 

31 

287 

1404 

30. 3N 

140.2E 

54^ 

Uppsala 

64 

329 

A39 

59. 5N 

U.llE 

16 

Vladivostok 

82 

358 

4577 

43. IN 

131. 5E 

55“ 

Vigna Di Valle 

44 

341 

6458 

42. IN 

12. IE 

>7 

Argentine island 

45 

319 

2310 

65. 2S 

64. 2W 

57 

Halley Bay 

40 

315 

1999 

75.35 

26. 4W 

19 

Bismark 

50 

349 

4946 

46. 5N 

100. 5W 

58 

Little America 

83 

318 

152 

78.0s 

162. OW 

20 

Caribou 

54 

370 

4098 

46. 5N 

68. OU 

62 

Port Aux Francais 

83 

375 

973 

49.2s 

70. 2E 

21* 

Edmonton 

55 

356 

6278 

53. 3N 

114. OW 

64 

Sterl ing 

42 

340 

1689 

38, 5N 

77. 3W 

22 

Green Bay 

49 

358 

4513 

44. 3N 

88. IW 

65* 

Toronto 

52 

362 

4298 

43. 4N 

79. iw 

23 

Moosonee 

57 

378 

1356 

51. 2N 

80. 4W 

66 

Ft. Collins 

39 

310 

I4l8 

40. 3N 

105. OW 

2k 

Resolute 

78 

390 

3878 

74 . 4n 

94. 5W 

67 

Boulder 

40 

332 

2577 

40.0N 

105. 2W 

26 

Aspendale 

36 

320 

6618 

38.0s 

145. OE 

68 

Belsk 

50 

341 

3971 

50. 5N 

20. 5E 

27* 

Brisbane 

24 

291 

5280 

27.3s 

153. OE 

69 

Hal 1ett 

41 

339 

4 00 

72.25 

170. IE 

28 

Dumont d'Urville 

85 

317 

333 

66 . 4 s 

140. OE 

70 

Mont Louis 

42 

336 

4450 

42. 3N 

2. IE 

29 

Macquarie tsl. 

50 

342 

4216 

54.3s 

158. 5E 

71 

Petoria 

15 

260 

1799 

25.5s 

28. IE 

30 

Marcus Island 

28 

270 

1970 

24. 2N 

153. 5E 

72 

Byrd 

45 

318 

947 

80. OS 

119. 3W 

31 

hauna Loa 

18 

276 

3372 

19. 3N 

155. 4W 

A 

73 

Ahmed abad 

16 

253 

3309 

23. ON 

72. 4E 

32 

Wei 1 ington 

43 

316 

1622 

41. 2S 

174. 5E 

74 

Varanasi 

18 

280 

4071 

25. 3N 

82. 5E 

34* 

Arhus 

63 

351 

5864 

56. IN 

10. IE 

75 

Dumdum 

17 

268 

2580 

22. 4N 

88. 3E 

35* 

Arosa 

45 

331 

4879 

46. 5N 

9.4E 

76 

Goose 

59 

380 

4752 

53. 2N 

60. 2W 

36 

Camborne 

49 

335 

2577 

50. IN 

5.2W 

77 

Churchi 1 1 

62 

387 

3572 

58. 5N 

94. OW 

38* 

Elmas/Cagl iari 

40 

331 

6011 

39. 2N 

9.OE 

79 

Ta] lahassee 

26 

306 

2432 

30. 3N 

84. 2W 



Table 10 continued 


Station 

Number 

Station Name 

0 

O3 

No. of 
Obser- 
vati ons 

Lati tude 

Longitude 

80 

Gan 

15 

264 

2215 

0 . 4 s 

73. IE 

81 

King Bedouin 

30 

330 

460 

70.3s 

24. 2E 

82 

Li sbon 

40 

301 

2235 

38. 5N 

9.IW 

84 

Darwin 

12 

264 

2706 

12.3s 

130. 5E 

85 

1 rkutsk 

93 

382 

3003 

52. 2N 

104. 2E 

86 

Karadag 

69 

302 

1927 

45. ON 

35. 2E 

87 

Kiev 

66 

338 

3880 

50. 2N 

30. 3E 

88 

Mi rny 

58 

315 

460 

66.3s 

93. OE 

90 

Ashkahabad 

58 

277 

4188 

37. 5N 

58. 2E 

91 

Buenos Aires 

28 

287 

1472 

34 . 4 s 

58. 3W 

92 

Hobart 

40 

327 

2241 

42. 5S 

147. 2E 

96 

Hradec Kralove 

51 

335 

3103 

50. IN 

15. 5E 

98 

Val Joyeux 

62 

304 

1801 

48. 5N 

2.0E 

99 

Hobenpei ssenberg 

46 

338 

1770 

47. 5N 

11. OE 

101 

Syowa 

47 

342 

892 

69.0s 

39. ^E 

102 

Bracknel 1 

^9 

352 

1512 

51. 3N 

0.5W 

103 

Albuquerque 

32 

297 

1493 

35. IN 

106. 4U 

104 

Bedford 

49 

357 

1569 

42. 3N 

71. 2W 

105 

Fai rbanks 

59 

383 

1288 

64. 5N 

147. 5W 

106 

Nashvi 1 le 

36 

33^ 

3836 

36. 2N 

86. 3W 

107 

Wallops Island 

38 

327 

1779 

37. 5N 

75. 3W 

no 

Huancayo 

9 

263 

3797 

12. OS 

75. 2W 

111 

Amundsen-Scott 

43 

325 

1377 

90.0s 

o.ow 

112 

Bolshaya Elan 

89 

364 

3118 

46. 5N 

142. 4E 

113 

Dushanbe 

64 

278 

3^54 

38 . 4n 

68. 5E 

115 

Kui byshev 

75 

330 

3328 

53. 2N 

50. 3E 

116 

Moscow 

79 

327 

2874 

55. 5N 

37. ^E 

117 

Murmansk 

87 

354 

2803 

68. 5N 

33. OE 

118 

Nagaevo 

102 

386 

2576 

59. 4N 

150. 5E 

119 

Odessa 

71 

329 

3430 

46. 3N 

30. 4E 


Station 

Number 

Station Name 

0 

O3 

No. of 
Obser- 
vations 

Latitude 

Long! tude 

120 

Omsk 

77 

369 

3'19‘t 

Slt-SN 

73. 2E 

121 

Riga 

74 

348 

3411 

56. 5N 

24. OE 

122 

Sverdlovsk 

69 

354 

4103 

56. 5N 

60. 4E 

123 

Yakutsk 

108 

366 

2668 

62. IN 

129. 5E 

128 

Ka raganda 

62 

269 

1019 

49. 5N 

73. IE 

129 

Pechora 

85 

303 

955 

65. IN 

57. IE 

130 

Petropavlovsk 

87 

357 

1277 

52. 5N 

158. 5E 

132 

Sofia 

41 

314 

921 

42. 5N 

23. 2E 

159 

Perth 

29 

295 

2314 

31.5s 

115. 5E 


Key Stations 
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Figure 48. Ozone Variations at 
Canada 
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Figure 50. Ozone Variations at 
Australia 
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Figure 49. Ozone Variations at 
Japan 
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Figure 51. Ozone Variations at 
India 
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Figure 52. Ozone Content, Each 
Station Weighted Equally 


Figure 53. Ozone Content, Each 
Hemisphere Weighted Equally 


Both data sets do show a decrease of ozone after '\^ 1970 and a 
distinct minimum in 1961. 

It is our conclusion that a carefully selected data set of 15 
stations indicates no obvious long-term trends in global total ozone. 
Because of the data sparsity over the oceanic regions and the strong 
bias toward the Northern Hemisphere (and especially Europe), we feel 
that analysts should utilize the total ozone data available with 
caution and careful inspection of parameters, such as the station a's. 

SATELLITE ANALYSIS 


Figure 54 indicates 100 days of total ozone data as measured 
by the Nimbus 3 IRIS sensor. These data have been latitudinal ly 
weighted to remove areal bias. The data extend through a period 
starting with the Northern Hemisphere spring (Southern Hemisphere 
fall) and ending with the Northern Hemisphere summer (Southern 
Hemisphere winter). The standard deviation for the data set is 
2.6 m atm-cm. 

During this period there was approximately 5% more total ozone 
observed by satellite in the Northern Hemisphere (318 m atm-cm) 
than in the Southern Hemisphere (303 m atm-cm) (Table 11). 
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Total ozone, m atm-cm 


Global average for period = 311.7 m atm-cm 
Standard deviation [o] = 2.6 m atm-cm 

Period deviation = 0.8% 


2 X 10® satellite analyzed data points 


15 20 25 30 
April 

5 10 15 20 

May 

25 30 5 10 15 20 25 30 

June 

Time — 

5 10 15 20 25 

July 


Figure 54. 

Total Ozone Data Measured by 
Nimbus 3 IRIS Sensor 







Table 11. Land and Sea Distribution of Total Ozone 

(April 16-Ju1y 22, 1969; Nimbus 3; 1.85 x lo 
data points; values in m atm-cm). 



Land 

Sea 

Total 

Standard 

Deviation 

Sea-Land 
(m atm-cm) 

Sea-Land 

% 

Global 

310.3 

312.1 

311.7 

2.6 

1.8 

0.58% 

No. Hemis. 

315.8 

319.1 

318.1 

7.8 

3.5 

1.04% 

So. Hemis. 

302.9 

303.1 

303.1 

8.9 

0.2 

0.07% 

30°N-60°N 

337.5 

344.1 

341.6 

14.3 

6.6 

1.96% 


Perry Gluckman presented his time series analysis which was 
done in the frequency domain rather than the time domain. His 
analysis tended to corroborate what other speakers had presented. 
Details of his presentation are not available for the proceedings. 

John Tukey of Princeton University discussed the use of 
exogenous variables. "While I do not for a moment undervalue 
physical insight or physical explanation, it is important to 
keep in mind that purely statistical considerations call for 
making adjustments of empirical size for any internally reliable 
exogenous variable that could possibly make sense. 




"We ought to do more to find and use exogenous variables." 

He concluded by suggesting some exogenous variables that might 
help: 

t Pittock's general circulation quantities 

• Reiter's energy sloshing and vacillation 

• Gluckman's intermonth adjustment to fixed dates 

t Gluckman's sector crossings - field reversal 

• ??? geomagnetic character figures 
Others suggested by conferees included: 

f local winds aloft 

• local barometer 

• local height of tropopause 

Tukey said, "Suppose we do adjust for the local barometer, then 
collect the global mean. Then we must think carefully about the 
interpretation if the mean barometric pattern is changing." London 
noted that this was already done, at least in part, in Pittock's 
general circulation. 

Tukey summarized his suggestion to "use the things we can 
trust--such as local pressure--and see what happens when we use 
them and then look for explanations." London responded that the 
key problem is the use of extra information in terms of filtering. 
"You are bringing up the key to the filtering problem in getting 
the real information." Tukey agreed saying we should "use all 
available principles of witchcraft and if some are roughly orthogonal 
we should use both." He restated Hill's methodology as using 
persistence and shocks to see what they tell us, then focusing the 
analysis on the shocks. This methodology, Tukey said, "does get 
you out of certain technical problems; it saves trouble with the 
data Hill had. If we can provide better data, perhaps he can do 
better." But Tukey concluded, "We cannot bypass Basher." 

James K. Angel 1, of NOAA, compared ozone trends with 
stratospheric water vapor, the temperature of the equatorial 
tropopause, and the north temperate latitude temperature (Fig. 55) 
to illustrate what he termed some "very interesting" results. 

Although the water vapor data record is short, beginning in 1964, 
and there are not many measurements (only one a month at the most), 
the total ozone is very well correlated with water vapor and 
distinctly out of phase with stratospheric temperatures. That 
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Figure 55. Comparison of Ozone Trends with Stratospheric Water 
Vapor, Temperature of Equatorial Tropopause, and North Temperate 

Latitude Temperature 
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is, the maximuin ozone occurs when the 
stratospheric temperature is the lowest, 
an unexpected result for which Reiter 
offered a meteorological explanation: 
stratospheric water vapor comes mostly 
from summer monsoons whereas ozone is 
a winter characteristic. He further 
suggested that the anticorrelation of 
temperature may be due to pressure 
distributions as suggested earlier by 
Pittock. 

Angel 1 expanded the puzzle with 
Umkehr data (Fig. S6) noting that, 

"If you accept Umkehr data, we see an 
increase instead of the expected de- 
crease due to CFMs. This problem is 
not really resolved but, where we 
should see a 5% decrease and we see 
instead 12% the other way, it makes 
us wonder." 
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Figure 56. Umkehr Data 


In subsequent discussion Angel 1 J 

pointed out that an anomaly in the data 
coincides with the eruption of Agung, 
leading him to question how Hill's 
analysis deals with such an anomaly. 

Hill said that he had misunderstood 
the previous question and that indeed 

the volcanic effect was in his analysis and that his techniques 
certainly try to quantify such interventions. Tukey elaborated 
that Hill's pre-whitening filter says nothing about mechanisms. 
"The effect of Agung is in there but it is not really significant 
— the whole question can be solved if Hill leaves out the quasi- 
biennial and uses only short analyses." 


Pittock agreed that using only short analyses would avoid 
"the primary problem of building in a prediction that Agung will 
happen again." 


Tukey, attempting again to summarize the issues at hand, said 
there are at least three parts to the problem — potential measurement 
troubles as described by Basher, the question of where you measure 
to gather global meaning, and the statistical factors— and each part 
"must be got at separately." One must allow for the real world 
because none of the other factors has yet included any natural trends. 


Referring to a three-dimensional chemical -dynamical model de- 
veloped by Jerry Mahlman of NOAA-GFDL, Campbell said this model 
showed poor agreement when used to compare the "true" global change 
in ozone to a global change estimated using Angell's best 53 sta- 
tions placing equal weighting on each of the 53 stations. However, 
Tukey interjected, "a very different idea" could be achieved by 
"sensibly weighting the stations geographically." Tukey suggested 
1) do a consistent (simplified) time-series analysis (short lags 
only) for say, 53 stations; and 2) study covariances and perhaps 
correlation coefficients between estimated shocks, and check the 
spectra, and some cross-spectra, of the estimated shocks. He 
recommended as further steps forward: "criticism" of empirical 

adjustment (regression) coefficients in terms of frequency bands 
(Fig. 57). 

Breakdown of such analyses as 



into at least a few frequency bands 


Figure 57. Empirical Adjustment Coefficient in Terms of Frequency 

Bands 

John DeLuisi asked if transformation can be made from historical 
ozone data to satellite data that, while global, will necessarily 
have some scatter. It would seem that a reasonable overlap would be 
at least one solar cycle. C. Desmond Walshaw noted that the Dobson 
instruments would be needed for some time. "Ozonesondes," he said, 
"were going to make the Dobsons obsolete and they did not." He 
continued, "Everyone who uses the total ozone network should be 
aware that there are all sorts of problems." (Fig. 58). 

The Dobson measurements are not only extremely important for 
the next 10 years but they are equally important as historical 
records if they can be corrected by "measurement archaeology." 

Komhyr noted that the basic problem in making Dobson spectro- 
photometer observations is the effect of pollution where it is 
estimated that errors of several percent can result. As far as 
NOAA's total O3 data are concerned, we have the basic calibration 
information that can be used to improve the quality of existing 
data; however, we do not have the necessary resources to make these 
corrections. 
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Figure 59. Worldwide Network 
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But after the current intercompan'sons Walshaw said, "We should 
have a satisfactory worldwide network for the first time." (Fig. 59) 

Joseph Drewry explained that his primary interest was mission 
analysis. "What future satellite missions do we need to determine 
total ozone?" he asked. 

"It is not obvious that time series analysis can address the 
ozone problem." Drewry proposed a possible solution of "letting 
the data develop a global spectral model of ozone in a natural 
coordinate system, and trying to minimize the variance of important 
model parameters with sampling analysis." Figure 60 shows the 
sampling capability of a simulated solar occultatlon mission over 
an ozone model based on weekly estimates over a 5“ x 15^^ global grid. 


13 WEEK MISSION DURATION (SPRING) 
H-600km. 1-50° 



Figure 60. Sampling Capability of Solar Occultation Mission 


He emphasized that the difference between the model estimates and 
estimates from the simulated mission reflects sampling distribution, 
not measurement errors. Drewry referenced Figure 61 when discussing 
empirical orthogonal functions as a technique for examining the 
information in a data set representative of global ozone data. 

He noted that "in the set of gridded data from which this example 
was taken, 98 percent of the variability about the monthly mean 
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can be explained by 6 eigenvectors; almost 85 percent is explained 
by the first principle component." He cautioned that data used in an 
analysis such as this could contain systematic errors which would 
be misinterpreted as ozone variability. 

Tukey suggested contouring the next six eigenvalues to get the 
last 2% and then "looking for the physics behind them, assuming there 
is some physics behind them." 

Glenn Brier presented a point of view as to how ozone helps 
understand the quasi -biennial oscillation using 26 years of data from 
Balboa (Fig. 62). He noted that, "If you look at a model with feed- 
back you should expect trends and, in a two-season system, you should 
get a biennial result." The actual result (Fig. 63) is very 
asymmetrical with respect to the seasons, yielding a picture of 
interventions and shocks which are not randomly distributed. 
Narasimhan Sundararaman mentioned that the high altitude pollution 
problem led his agency to modify the original question, "Can we find 
a trend in the ozone?" into a new, more-to-the-point question, "What 
is the optimum Dobson network that would really give us the trends 
and how do we get that network?" 
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To get a deeper understanding of the measurement error problems 
for ozonesondes and Dobson instruments, Heath suggested developing 
independent measurement checks, and comparing high quality Dobson 
data with satellite information. Rocket measurements, too, are 
possible. Although quite difficult, a standard rocket payload has 
now been developed. The rocket program, begun 10 years ago, should 
provide information on the total ozone trend by noting a trend in 
the 40-60 km altitude range. 
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SESSION IV: CONCLUDING REMARKS 


At the conclusion of the Symposium, the chairmen from Sessions I and II 
were given an opportunity to make summary statements. 

Dr. Hill began by thanking the conveners of the Symposium for the 
opportunity to present details of his time series analysis and the fruitful 
interchange that resulted. He reaffirmed his belief in empirical methods 
as "letting the data speak for themselves" rather than interjecting into 
models preconceived physical mechanisms that may not be supported by the 
empirical evidence. He conceded that empirical methods can lead to physically 
meaningless or unexplained results and, therefore, must be interpreted in 
light of plausible physical mechanisms. Dr. Hill concluded by expressing his 
hope that the dialogue begun at this symposium will continue. 

In his summary remarks. Dr. London stated that there seem to be no serious 
objections to the statistical methods used. It is only the conclusions that 
are questioned, on the grounds that (1) the length of record was probably too 
short to eliminate low frequency effects of meteorological variabilities, 

(2) there may be systematic long-term trends affecting the observational 
system (giving incorrect data variations), and (3) stations chosen for the 
trend may not be representative of their geographic area and, therefore, would 
not give a correct global average. 

Summarizing the suggestions offered during Session III, Dr. London under- 
scored the recommendation that the same statistical methods be applied to 
meteorological data for which there are long, compatible series (e.g., tempera- 
ture, precipitation, drought index, etc.) and where known trend changes have 
taken place (e.g., change from Northern Hemisphere warming to cooling around 
1940). A second suggestion was that further research and data "washing" be 
done to make the various observational series homogeneous. The effects of 
optical wedge deterioration, atmospheric aerosol variation, solar irradiance 
variation, etc., need to be evaluated with more precision than has been done so 
far. "It should be emphasized that the importance of the problem dictates 
that reasonable sums of money must be expended to support this type of research." 
Finally, referring to the geographic representativeness of the data. Dr. London 
emphasized that a coupled satellite ground-based observational system is 
required to determine global long term trends. This requires maintenance and 
improvement of the Dobson network and long term planning for a satellite 
observing system. 

Dr. London concluded by thanking the NASA sponsors, in particular 
Dr. Greenwood, for convening the Symposium, and the attendees who took time 
from their busy schedules to participate. 

Dr. Greenwood also thanked the attendees and suggested that the partici- 
pants send him their comments and/or recommendations after they have had time to 
reflect on the discussions. "A role that NASA can play is to encourage a 
continuing dialogue and we are open to suggestions on how best to do this." 
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